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I. INTRODUCTION 

The influence exerted by zero point ionic motion on the 
stability of the Peierls and spin-Peierls lattice distorted 
states enters as a key ingredient in the elaboration of a 
general theoretical description of these phases. Quan- 
tum fluctuations are known to cause a downward renor- 
malization of the order parameter and the corresponding 
electronic gap, if not their complete suppression as it 
is the case for spin-Peierls order. One is confronted to 
such situations in low dimensional conductors and insu- 
lators for which the characteristic phonon energy is not 
only finite in practice, but may exceed by far the tem- 
perature scale at which the lattice instability takes place. 
These cases are exemplified in spin-Peierls systems like 
the inorganic compound CuGe03fi^ the organic system 
MEM(TCNQ) 2) 2 and also members of the (TMTTF) 2 X 
series of organic compounds for which non adiabaticity 
emerges as one moves along the pressure scale, giving rise 
to quantum criticality for the spin-Peierls transition. 3 

The first systematic studies of quantum effects on the 
Peierls-type distorted states go back in the eighties with 
the world-line Monte Carlo simulations of Hirsch and 
Fradkin4^ These simulations were made on the one- 
dimensional tight binding and Holstcin electron-phonon 
models, also known as the Su-Schricffer-Heeger (SSH) 6 
and the molecular crystal (MC)£ models. The stability 
of lattice distorted phases was determined as a function 
of the ionic mass and the strength of electron-phonon 
coupling. The phase diagrams of the models were out- 
lined for both spinless and spin-i fermions at half-filling. 
These initial works were followed by a variety of numeri- 
cal techniques applied to the same models and extended 
to include direct interactions between fermions. That is 
how density matrix renormalization group (DMRG)r^i2i 
exact diagonalizationsjii and quantum Monte Carloi^ 
techniques to mention a few, have contributed to provide 
a fairly coherent picture of the influence wielded by zero 
point lattice fluctuations in one-dimensional electron- 
phonon systems. 

On the analytical side, these progress were 



preceded 1 ^ and accompaniedi^^^^^^^ by 
a whole host of approaches applied to study retar- 
dation effects on lattice distortion at intermediate 
phonon frequencies. The renormalization group (RG) 
metho d 15 ' 17 i 18 ' 21 i 22 has been one of the routes proposed 
to deal with this problem. A variant of the RG method 
will be further developed in this work. Our analysis 
starts with the effective fermionic formulation of the 
electron-phonon problem, which is expressed in terms 
of a fermion gas in the continuum with weak retarded 
interactions. Such a formulation for the SSH and MC 
models has been investigated long ago by the two-cutoff 
scaling method. 15 In this approach the characteristic 
bandwidth energy Eq for fermions and the vibrational 
energy lu c (H = 1) for phonons determine the form of 
flow equations for the electronic scattering amplitudes, 14 
whose singularities signal the creation of gaps and 
long-range order at half-filling. Thus when the electronic 
mean- field energy gap Ao — emerging below Eq in the 
adiabatic weak coupling theory — is larger than uj c , 
quantum corrections are neglected and the flow is equiv- 
alent to a ladder diagrammatic summation compatible 
with the unrenormalized static scale Ao for the gap. On 
the other hand when Ao < oj c , the scattering amplitudes, 
though still governed by the ladder flow down to lu c , 
are considered as effective unretarded interactions at 
lower energies. Below oj c the flow becomes impregnated 
by vertex corrections and interference between different 
scattering channels. In accord with the well known re- 
sults of the one-dimensional electron gas model ) 23 ! 24 ' 25 ! 26 
the classical gap Ao is then an irrelevant scale and the 
system enters in the non adiabatic quantum domain 
where either a gapless or an ordered massive phase can 
occur. 

While the two-cutoff RG analysis can provide a sim- 
ple and reliable criteria to map out the essentials of 
the quantum-classical boundaries of the phase diagram 
for both models in the weak coupling sector^ it says 
nothing, on the other hand, on how the gap varies over 
the whole phonon frequency range. This is not only 
of practical importance, when e.g., the theory is con- 
fronted to experiment in concrete cases, but also clearly 
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needed on general grounds when one raises the ques- 
tion of the nature of quantum-classical transition as a 
function of phonon frequency. This drawback is not a 
weakness of the RG method in general but rather en- 
sues from the frequency dependence of couplings, which 
in the two cut-off scaling approach, barely reduces to 
the minimum found in either the adiabatic or non adia- 
batic limit. A continuous description of retardation ef- 
fects would require that the full functional dependence 
of scattering amplitudes on the frequencies be restored, 
a possibility that can be liken to what has been done 
in two-dimensional and quasi-onc-dimcnsional RG for 
the functional dependence on scattering amplitudes on 
the momentum . 27 ! 28 ' 29 ' 30 ' 31 ! 32 Very recent progress along 
these lines show that it is indeed a promising avenue.— 
In this paper we shall revert to the RG approach as 
developed in Refs . 26 ' 34 and extend its formulation to in- 
clude the frequency dependence of scattering amplitudes 
introduced by the electron-phonon interaction. We re- 
visit the classical and quantum aspects of fermion driven 
lattice instabilities. Our analysis is done at the one-loop 
level and covers the gap determination and the structure 
of the phase diagram of the MC and SSH models for 
both spinless and spin-i fermions. Although the gener- 
alization to incommensurate band filling and situations 
where the direct Coulomb interaction is included would 
cause no difficulty, we have restricted our analysis to re- 
tarded interactions at half-filling. In Sec. II, we introduce 
the electron-phonon models and recall the derivation of 
their respective bare retarded interactions in the frame- 
work of an effective fermion gas model. We pay special 
attention to the SSH model in the spinless case in or- 
der to include the momentum dependent umklapp term 
to the interaction parameter space, which is so impor- 
tant for long-range order of this model. In Sec. Ill the 
one-loop level flow equations for the retarded scattering 
amplitudes and response functions are derived for spin- 
less and spin-i fermions. In Sec. IV we compute the 
variations of the gap over the whole frequency range and 
discuss the structure of the phase diagram and the na- 
ture of the quantum-classical transitions for the MC and 
SSH models. We conclude in Sec. V. 



II. THE MODELS AND THE PARTITION 
FUNCTION 

A. Models 

The one-dimensional electron-phonon models that we 
shall study using the RG method are the MC and the 



SSH models. The MC model describes the coupling of 
fermions to optical molecular phonon modes, whereas 
for the SSH model the electron-phonon interaction re- 
sults from the modulation of electronic energy by acous- 
tic phonons. In Fourier space, the two one-dimensional 
models Hamiltonians can be written in following form 

H = H n + H p h + H i 

= E e ( fc ) c Lc fe , CT +^^(^ + -) 

k,a q 

+ L ~ h l)A+ q ,* c kAbl + 6_,). (1) 

Here Hq is the free fermion part and e(fc) = — 2tcosk is 
the tight-binding energy spectrum with t as the hopping 
integral (the lattice constant a = 1 and L is the num- 
ber of sites). c\ (ck,cr) creates (annihilates) a fermion 
of wave vector k and spin a. H p h and Hi terms corre- 
spond to the free phonon and electron-phonon interaction 
parts, respectively, and in which b* q (b q ) creates (annihi- 
lates) a phonon of wave vector q. For the MC model, 7 
the intramolecular phonon energy and the interaction are 
given by 

UJ q = UJ , (2) 

g(k,q) = A /v / 2Mo^, (3) 

which arc both independent of the momentum. Here 
Ao > is the amplitude of the electron-phonon in- 
teraction on each molecular site whereas the frequency 
lo = ^/k /M is expressed in terms of the elastic con- 
stant Ko and the molecular mass Mq. 

For the SSH model,— the corresponding quantities read 

uj q = cj D | sin - 1 , (4) 

q(k,q) — i4 ■ ^ D sin f cos (k + , (5) 

y ' x/2M n u n 2 V 2 h W 



where lob — 2\J kd /Md is the acoustic phonon energy 
at q = 2k p, namely at twice the Fermi wave vector kp — 
ir/2 at half-filling. is the ionic mass and kd is the 
constant force of the one-dimensional lattice. 



B. The partition function 

Following the trace over harmonic phonon degrees of 
freedom in the interaction Matsubara time representation 
of the grand canonical partition function Z, one can write 



Z = Tr c e-P""-^ Tr ph e-^ h T T exp j -jf" Hi(r)dr j 
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Z ph Tr c e -PBo-»N Tr exp | _ 1 f f 9(k, q)g(k', -q)D(q, r - r')c{ +q 



-q,a 



,(r')c kl ^(T')c kt<T (r)dT'dT\ 



(6) 



where Z p h is the partition function of bare phonons. 
The phonon integration introduces an effective 'retarded' 
fcrmion interaction mediated by phonons and described 
by the bare propagator 

D(q 7 T-T r ) = e -^- r 'l +2(e? u * -ly 1 cosh(w g (t-t')) . 

The remaining trace over fermion degrees of freedom can 
be recast into a functional integral form 

Z = Z ph JJ®iP*®iP e s[r ^\ 

= z ph JJn^Dip e so[rM+si[r-M^ ( 7 ) 

over the anticommuting Grassman fields ip. In the 
Fourier Matsubara space, the free fermionic action is 

Soir, = £ KWY^A^Ah ( 8 ) 

p,k,<7 



where 



G° p (k) = [iu - e p (k)] 1 (9) 



is the bare fermion propagator for 
k = (k, u) = ±7rT, ±3irT, . . .) (ks = 1)- The fermion 
spectrum e(k) — /i sa e p {k) — v F (pk — kp) is linearized 
around the right (p = +) and left (p = — ) Fermi 
points ±kp. The bandwidth cut-off Eq = 2Ep 
is twice the Fermi energy Ep — vpkp. The inte- 
gration of the fermion degrees of freedom becomes 

The interacting part 57 of the action reads 



T 
2L 



{p,fc,cr} 



+2,1^3+4! 



(10) 



where momentum conservation is assured modulo the re- 
ciprocal lattice vector G = ±4kp, allowing for umklapp 
scattering at half-filling. In the Fourier-Matsubara space, 
the interaction takes the form 

g{ki,k 2 ;k 3 , h) = g(ki,k 3 - k 1 )g(k 2 , k 4 - k 2 )D(k 3 - fei), 

(11) 

where 

D(h-k) = -2 - 2 ^ , 

is the bare phonon propagator. We can now proceed to 
the 'g-ology' decomposition of this interaction. This will 
be done separately for fermions with and without spins. 

In the first place, for spin-i fermions, we shall con- 
sider the three standard couplings between fermions on 
opposite Fermi points 

3l(wi,CJ 2 ,W 3 ) EE g(±k F ,0Ji,=fk F ,U 2 ;TkF,U3,±kF,U>4,), 

#2(^1,^2,^3) = g(±k F ,u}i,Tk F ,uJ 2 ;±kF,cj 3 ,TkF,V4), 
cj 3 (ll>i,uj 2 ,u} 3 ) = g(±k F ,0Ji,±k F ,u> 2 ;Tk F ,u 3 ,TkF,U4), 



(12) 

for retarded backward, forward and umklapp scattering 
amplitudes, respectively (here the forward scattering of 
fermions on the same branch is neglected) . According to 
Eq. ([3]), the bare frequency dependent couplings for the 
MC model become 

54=1,2,3(^1,^2,^3) = — — f 2 , (13) 

where 5^=1,2,3 = — Aq/ko is the (Mo-independent) attrac- 
tive amplitude. Similarly for the SSH model, one has 
from Eq. © 

gi 3 

31,3(^1,^2,^3) = — — 2 2 , (14) 
1 + w 3 _ 1 /o; £ , 

where the amplitudes gi, 3 = z f4A|,/k£) are also Mp>- 
independent. For the SSH model, the bare forward scat- 
tering amplitude g 2 vanishes for the exchange of zero mo- 
mentum phonon but it will be generated at lower energy 
by the renormalization group transformation. 

For spinless fermions, the backward scattering is in- 
distinguishable by exchange from the forward scattering 
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and both can be combined to define an effective forward 
scattering term of the form 

- ,g(±fc F , u 2 , T&f , ; =f^f , ^3 , ±&f , w 4 ) 
52 51 (15) 



1 



-l/ W 0,Z 



1 



J 3-2 



1^1,1 



where the mass independent amplitudes are 
gi,2 = — Aq/kq for the MC model, and g\ — —A\ 2 d /k,d 
and gi = in the SSH case. 

As for the umklapp scattering in the spinless case, it 
must be antisymmetrized with its own exchange term to 
give the following two contributions 



1 



^[gih, k 2 ; k 3 , k 4 )-g(k 2 , h; k 3 , k 4 )} k k lZ k k 2 = 53(^1,^2,^3) 



+ ,9m(wi, u>2, W3)(sin k\ — sin fc2)(sin k 3 — sin k 4 ). (16) 

The first contribution corresponds to a local umklapp 
term defined for incoming and outgoing fermions at the 
Fermi points. It takes the form 

.93(^1,^2,^3) = ^ [g(±k F ,u} 1 ,±k F ,uj2- 1 TkF,U3,TkF,uJ4) 

- g(±k F , u 2 , ±k F , u\\ ^fk F ,u>3, T-^f, Ui)] 
= 93 (17) 

1 + Ui^/ul D 1 + W3-2/ w 0,d' 

This term is present for both models, where 173 = — Aq/ko 
is attractive for the MC model, and 53 = 4A|,/k_d is re- 
pulsive for the SSH model. The second term of (TT6"]) is a 
non local - momentum dependent - umklapp contribu- 
tion and is only present for the SSH model. Actually this 
additional contribution follows from the antisymmetriza- 
tion of Eq. (fTTjl and the use of ([5|) under the permuta- 
tion of incoming and outgoing frequencies and momen- 
tum (these last, not at the Fermi points). Its frequency 
dependent part reads 



g u (oji, uj 2 , L03) = g u 



1 



1 



tot 



- U) 



3-2 



«2(4) 



(18) 

\ 2 D /2n D . From 
the leading k- 



where the amplitude is given by g v 
Eq. (TT6")) . it follows that for k^ <~ 
dependence of the non local umklapp is cx (fci — ^2X^3 — 
ki) which has a scaling dimension of -2. This term is 
therefore strongly irrelevant at the tree level but becomes 
relevant beyond some threshold in the elect ron-phonon 
interaction. Such umklapp contributions are well known 
to play a key role in the existence of long-range order for 
interacting spinless fermions ) 35 ' 36 as it will show to be 
the case for the SSH modelJ.&±£ 



III. THE RENORMALIZATION GROUP 
TRANSFORMATION 

The renormalization group transformation for the par- 
tition function will follow the one given in Ref i 26 i 34 One 



then proceeds for Z to the successive partial integration 
of fermion degrees of freedom, denoted by $W having 
the momentum located in the outer energy shells (o.s) 
±Eo(£)d£/2 above and below the Fermi points for each 
fermion branch p. The remaining (<) degrees of free- 
dom are kept fixed. Here Eq(£) = Eoe~ e is the scaled 
bandwidth at £ > 0. The integration proceeds by first 
splitting the action S — > S[ip*, + Sq + Si into an in- 
ner shell part at £ and the -tp - dependent outer shell 
terms Sq and Si. Considering Si as a perturbation with 
respect to the free outer-shell action Sq, the partial inte- 
gration at the one-loop level is of the form 



S[ip',ip]e 



2)^*2)^ e S[^k+(Si>o. 8 +i<s?>o. s +....( 19 ) 



Here the interacting part is made up of three pertinent 
terms, i.e., Si = Sf 2 + Sf 2 + Sj 2 , for all possibilities 
of putting simultaneously two outer shell fields in the 
2k F electron-hole Peierls channel (Sf 2 ~ '4 ! X4 ! -i>~'^+ + 
V'+V'+V'-V'- + ■■■), the zero momentum fermion-fermion 
Cooper channel {Sf 2 ~ ip*^ip*lip-4>+ + ■••), an d the Lan- 
dau channel (Sf 2 ~ + ■■■)■ 

The lowest order outer shell statistical average ( Si ) . s 
comes from the Landau part and gives rise to the self- 
energy corrections 5Y.(lu) of the one-particle Green func- 
tion, which becomes G" 1 + z<5£(cj). As for the contrac- 
tions \ ( Sj )o. s , only the singular Peierls and Cooper scat- 
tering channels are retained; with four fields in the inner 
shell, these correspond to corrections to the coupling con- 
stants. Both corrections define the renormalized action 
S[ip*, ip]e+M at the step £ + d£. 

The evaluation of outer shell contractions (Si) . s for 
the self-energy at the one-loop level leads to 

(Si)o.s = iST,(u) ^2 ^p,a(k)ip P ,a(k), 



<5£(w) = -kv f -^2 /2 9s(u',u,u)G-(k,u'). 



w' {fe}o 



(20) 



In the low temperature limit, the flow equation for the 
self-energy becomes 



f +ca du' r 

/ IT- \9s{u ,u,uj) 
J-00 27r 1 

„ (E (£)/2)(u,'-Z(u,')) 



(u/ - £ p (u/)) 2 + (£oW/2) : 



where 



(21) 



(22) 
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for spin-i fermions, and 

g s (u',uj,uj) = -g f (uj,uj',uj), (23) 

in the spinless case. Here g~i({u>}) = gi({uj})/nvp, di 
ill. and E(w) = at I = 0. 



[3((^ 2 ) 2 >o. s .] and Cooper [i((S' / c 2 ) 2 )o.s.] channels. For 
the MC and SSH models defined by the initial couplings 
(|13H14j) . these interfering contractions lead to the flow 
equations 



A. RG flow for couplings: spin-| fermions 

The one-loop contractions ^(Sj) . s amount to eval- 
uate the outer shell contributions of the Peierls 



d t g x (uj x ,uj 2 ,uj 3 ) = 



dw ( 
2ttI 



dig 2 (uj u uj 2 ,uj 3 ) = / — j 



; + °° dw r 

dig 3 (uj x ,uj 2 ,uj 3 ) =2 \ — | 



2 5i (w, uj 2 , w + w 3 - uj x )g x (uj x , uj + uj 3 - w x , uj 3 )I p (w, w 3 - wi) 

§ X (uj, W 2 , Cl> + W 3 - Wx)5a(w + w 3 - ^1,^1, 103)lp(b), uj 3 - UJ X ) 
g X (uJ X ,UJ,UJ 3 )g 2 (uJ 2 ,UJ + UJ X - UJ 3 ,Uj)I P (uJ,UJ X - UJ 3 ) 

2g 3 (uj 2 ,uj,uj x +uj 2 - u) 3 )g 3 (u>i, uj + uj 3 - uj x ,uj 3 )I p (uj,uj 3 - uj x ) 
g 3 (uj x ,uj,u> 3 )g 3 (uj + uj x — uj 3 ,uj 2 ,uj x + uj 2 - uj 3 )I p (uj, uj x - uj 3 ) 
g 3 (uj 2 ,uj,uj x + uj 2 - u> 3 )g 3 (u> + uj 3 - uj x ,uj x ,uj 3 )I p (uj,uj 3 - uj x ) 
g 2 (uj 2 ,uj x ,uj + uj x + uj 2 )gi{-uj,uj 2 ,uji + ui 2 + uj) 

g~\{u)l, UJ 2 , UJ + LOl + L0 2 )g 2 (uj + UJ X + LU 2 ,—LO, L0 3 ) Ic(w, UJ X + W 2 ) j 

g 2 (u>i , uj + u> 2 - u> 3 , uj)g 2 (uj , u> 2 , u> 3 ) 

g 3 {uji,u) + uj 2 - uj 3 , uj)g 3 (uj, u> 2 , uj 3 ) Ip(uJ, ui 2 - uj 3 ) 

g 2 (uj x ,uj 2 ,uj + uj x + uj 2 )g 2 {uj + uj x + ui 2 , -uj,uj 3 ) 
g x (uj 2 , wi, uj + UJ\ + UJ 2 )g x (—uj,uj + uj x + uj 2 ,ui 3 ) Ic(uj,uJi + UJ 2 )\ 

2g x (uj x , uj + uj 3 - ujx,uj 3 )g 3 (uj,uj 2 ,uj + uj 3 - uj x )I p (uj,uj 3 - uj x ) 
g x (uj x ,uj,u; 3 )g 3 (uj + uj x -uj 3 ,uj 2 ,uj 1 + ui 2 - uj 3 )I p (uj,uji - uj 3 ) 

g 3 (uj x ,uj,uj 3 )g 2 (uj + uj x — uj 3 ,uj 2 ,uj x + uj 2 - uj 3 )I p (uj,uj x - uj 3 ) 

g 3 (uj,ui 2 , uj 3 )g 2 (ui + ui 2 - uj 3 ,uj x ,uj x +ui 2 - uj 3 )I p (ui,ui 2 - uj 3 )\. 



(24) 



(25) 



(26) 



The momentum shell Peierls and Cooper loops I P (uj,tt) frequency f2 are by using the Green function with self- 
and Ic(uJ,&) at internal uj and their respective external energy corrections 
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I P {u,n)d£ = — T~ G + {k + 2k F ,uj + n)G^{k,uj), 

{k} 3 

d£ " (^-EH)(^ + r>-E(^ + r>)) + iff 2 (i) 

2 oU [( W - E( w ))a + ±E 2 (*)][(u; + - E(w + n))» + ±E 2 (£)] ' 
I c {u,VL)dl = ^ ^ G + (fc,w + fi)G_(-fc,-a;), 

{fe}o s 

= ^ m - sM)(" - n + s(n - ^)) + |£ 2 (i) 

2 01 j [( W - EH) 2 + \E 2 {i)][(u - n + £(fi - c)) 2 + l£ ^)] ' 



B. RG flow for couplings: spinless fermions 

Owing to the nature of umklapp scattering which in 
the spinless case is different for the MC and SSH mod- 



I 

els, we shall proceed separately for each model. Thus 
for the MC model with a local umklapp term, the outer 
shell contractions \{{Sf 2 ) 2 )a.s. and ^((Sf 2 ) 2 } .s. for the 
Peierls and Cooper channels allow us to write 



+°° doj 
2 



-{ gt(wx,(jiJ + (jJ2-U3,u)9f(v,V2,va) 

TT v. L 

+ 53(^1,^ + ^2-^3,^)53(^,^2,^3) Ip(u,W2- OJ3) 

- (jf(uJl,UJ2,LU + LUl + LU 2 )gf(0J + U\ + U>2, —u},u}3)Ic(u,ui + w 2 )|, 



^33(^1,^2,^3) = 2 J ^ I g 3 (uji,u,uj3)gf(u + u>x - u} 3 ,u 2 ,uji +u 2 - w 3 )Ip(cj,Wi - u> 3 ) 

+ 33 (^->, ^2, ^3)3/ + ^2 - W3,0;i,Wi + CJ 2 - W 3 )/p(u;,£J2 ~ W 3 )j, 



(29) 



(30) 



which are subjected to the initial conditions (fT"5f and (fT7|) 
of the MC model. 

The flow equations for the SSH model presents 
an important difference because of the additional 
k— dependent umklapp term (fig]) . At variance with 
gf and 33, this coupling acquires a non zero scaling 
dimension at the tree level. Therefore the momen- 
tum, energies and fields must be rescaled after each 
partial trace operation in Eq. (|19|) . which restores 
the original bandwidth cutoff. Thus following the 



outer shell integration, one applies the transformations 
k' = sk, u)' = suj, ipW = s" 1 / 2 ^*), T i = sT; jj = s -i Lj 
and oj' d — sujd (M' d = s~ 2 Md, the spring constant Kd 
is kept fixed), where s = e dl . This gives the scaling 
transformations for the outer shell corrected couplings, 
namely g'f 3 ({oj'}) — s°gf, 3 ({uj}) for the local part and 
g' u {{uj'}) = s~ 2 g u ({uj}) for the non local part. The flow 
equations for the SSH model then become 



digf(oJi,Lu 2l uj 3 ) = 



+ °° dco 



duo r 

27r\ 



gf{w\, LU + U2 - ^3,^)3/(^,^2,^3) + 53(^1, cj + u> 2 -0^3,^)33(^,^2,^3) 

+ g u (tJx,U + (J2 - ^3,^)33(^,^2,^3) +33(^1,^ + ^2 - UJ 3 ,L0)g u (uj,UJ2,C03) 
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+ g u (wi,w + w 2 - u 3 ,u)g u (uj,uj 2 ,uj 3 ) I p (u,lj 2 - oj 3 ) 

- gf(u)i,ui2,u + u)i + u)2)gf{u + ui +W2,—u),(jJa)Ic(v,vx + u;2)\, (31) 



dtg 3 {u\, lu 2 , W3) = 2 



+ g 3 (uj,uj 2 ,uj 3 )gf(uj + uj 2 - W3,u>i,ui + lu 2 - lu 3 )I p (uj,uj 2 - W3)}, 



(32) 



?u(wi, w 2 , u> 3 ) 



2g u (uj 1 ,uj 2 ,uj 3 ) 
50 f 



LU,UJ 3 )gf(uJ +U)l— UJ 3 ,UJ 2 ,U!l +uj 2 - UJ 3 )Ip(uJ,UJl - UJ 3 ) 

+ g u (uj,uj 2l uj 3 )gf(uj + uj 2 - W3, wi, Wi + W2 - uJ 3 )Ip{bj,bJ2 - W3)}, 



(33) 



which are subjected to the initial conditions (fT5|) , (fT7|) . 
and |g5D . 



C. Response functions 

Staggered 2/c^ density- wave and zero pair momen- 
tum superconducting susceptibilities can be computed by 
adding a set of linear couplings to composite fields {h^} 
in the bare action at £ = 0. This gives the source field 
term 

s h [p,,i>] = £ [ E <(nX(w,w + n)o£>,fi) 



+ c.c 



(34) 



where zf^ and the renormalization factors of 



'i'p 



the corresponding source fields, with zjxfX = 1 for the 
boundary conditions at £ = 0. For spin-i fermions, we 
shall focus on the 2kp susceptibilities for 'site' M = + 
and 'bond' M = — charge (CDW, BOW: fi p = 0), and 
spin (SDWaj^z, BSDW^j^z: ^=1,2,3) density-wave cor- 
relation of the Peierls channel. The corresponding com- 
posite fields are 



0£>,fi) = I [0^(0;, Sl) + MO; p (u, -O)], (35) 



Efe a/3 - 2kp,UJ - n) 



(36) 



In the Cooper channel, we consider the uniform super- 
conducting singlet (SS: /i c = 0) and triplet (TS Xj y iZ : 
/i c = 1,2,3) susceptibilities. The corresponding com- 



posite fields at zero pair momentum are given by 



k,a/3 



(37) 

are the Pauli 



For both channels, <jq — 1 and 012,3 = 
matrices. 

In the case of spinless fermions, only the 2kp 'site' 
CDW and BOW susceptibilities survive with 



1 T 



^2[ip_{k - 2k F ,u> - n)i>i(k,u) 
k 

+ Mipl(k + 2k F ,u-Sl)ipt(k,oj))]. (38) 

In the superconducting channel, only one susceptibility 
is considered with the corresponding pair field 



Me (u,Sl) = Jy5^V-(-*.-w+n)^+(A,w). (39) 



Adding (|34|) to the action S in ([7]), the partial inte- 
gration (TIT)]) at the one-loop level yields additional outer 
shell contributions that correct Sh and which gives the 
recursion relation 

Sh[ip*,ip]i+<u = S h [tp*,tp]i + (S h Si) . s + -(Sl)o, s + ■■■ 

(40) 

The second term (ShSi) . s is proportional to O^/i* and 
its complex conjugate, and leads to the flow equations for 
the renormalization factors z„ of the pair vertex parts. 
In the density-wave channel, its evaluation leads to 



x [g„ p {u,u\u> + ri)+Mg 3 (cj,uj' + n,u')]l P (cj',tt)}, 
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where 

5 Mp =o(^, w', w + si) = 5 2 (w', w, w + n), 

-2gi(w, + 0), 

.g Mp #o(^, w', w + Q) = g 2 (u', uj,oj + SI), 

Q 3 (uj,uj' + Sl,u}') = g 3 (id,id' + SI,uj') 

— 2g 3 (ui, u' + SI, u -\ 

for fermions with spins, and 



SI), (42) 



(41)advantage of certain symmetries namely, the time inver- 
sion, left-right Fermi points symmetry, and the exchange 
symmetry between the incoming (wi,o>2) and outgoing 
lu^. == uj\ -\- lU'2 — ^3) frequencies. The last symme- 
try antisymmetrizes the initial conditions for the spin- 
less fermions case, especially for the umklapp process g 3 . 
We thus have to calculate 932 different functions for each 
gi. The same procedure is used to calculate the response 
functions and susceptibilities. The flow equations are nu- 
merically solved until the most singular susceptibility di- 
verges with the slope irvpdexff = 10 6 , which determines 
the critical value £ c at which the algorithm is stopped. 



3mp( w > w '' w + SI) = g f (w,w',w + Sl), 

Q 3 (uj,uj' + SI,uj') = —g 3 {u),u)' + Sl,uj'), (43) 

in the spinless case. Similarly, for the superconducting 
channel, one gets 



d e z^ c (uj,~-uj + S2) = 



( °°^L{ Zflc (uj',-uj' + Sl) 
x g Me (u, (J, Sl)I c {u', SI)}, (44) 



where 

g^ c= o{uj,uj',Sl) 



—gi(Sl — u, u>, u r ) — 92(Wj SI — u>, u>'), 

gx(Sl-oj,oj,oj') -g 2 (uJ,Sl~ uj,uj'). 

(45) 



for spin-i fermions, and 



gn c (u,u',Sl) = -gf(w,Sl-u,u') 



(46) 



in the spinless case. 

As a result of the partial trace integration, the last 
term of (|40[> . which is proportional to h^lth^fls , is gen- 
erated along the flow and corresponds to the susceptibil- 
ity in each channel considered, namely 



<!,(«)= 0"* )- 1 /_ +0 °S{l«,(^- + «)l 



x(2s + l)/ c(P) (o;,Q)}, 



(47) 



which has been defined positive (x^fl-, (SI) = at £ = 0), 
and where s is the spin. 

We close this section by a digression on the numerical 
aspects associated with the solution of the above equa- 
tions. Their numerical evaluation makes use of patches 
in the frequency manifold. The frequency axis is dis- 
cretized into a total of 15 subdivisions or patches be- 
tween the maximum values w max = ±1.5 Ep, which serve 
as bounds of integration for the frequency. The interac- 
tion is taken as constant over each patch where the loop 
integrals are done exactly. In order to reduce the num- 
ber of frequency dependent coupling constants, we take 



IV. RESULTS 
A. Adiabatic limit 

The results at non zero-phonon frequency will be com- 
pared to those of the adiabatic limit where i^o(D) —> 0. In 
this limit, the initial conditions given in Sec. [II] for both 
models show that either ^3-1 — > or o->3_2 — * 0, indicat- 
ing that no phonon exchange between fermions at finite 
frequency is possible. In the spinless case the gf coupling 
Eq. (fT"5|) reduces to its backward scattering part. There- 
fore only close loops contribute to the renormalization 
of both the coupling constants, susceptibilities, and one- 
particle self-energy S; the latter being vanishingly small 
in the adiabatic limit. 

The flow equations for fermions with s = 1/2 (resp. s — 
0) (|24j|26j) [resp. Eqs. (|29H30j) ] can be recast into equa- 
tions for g\ and 173, which become independent of fre- 
quencies 



dt{g\ ±53) 



-{2s- 



l)(~ gi ±~g 3 ) 2 /2. 



The solution is obtained at once 



ffi(*)±fls(i) = 



91 ±53 



l + H2 S + l)(~ gi ±g 3 )f 



(48) 



(49) 



which presents a singularity at £q = — 2[(2s + l)(gi ± 
g 3 )] _1 for combinations of bare attractive couplings gi ± 
g 3 found in the MC (+) and SSH (-) models (Eqs. 
(p ]) -rT4| ) . and (j!5til6|) V This signals an instability of the 
fcrmion system and the formation of a Peierls state with 
a — mean-field (MF) — gap Ao(= 2Epe~ l °), which takes 
the BCS formic 



A = 2Ep exp 



(-2/(25 + 1)1^ 



(50) 



This singularity is present in the pair vertex factors z^ _ (l 
at SI = in either CDW or BOW channel depending of 
the model. In the adiabatic limit this can be seen by 
retaining only closed loops in (|41[) . where for frequency 
independent couplings, zjf _ becomes in turn indepen- 
dent of u> and obeys the following flow equation at SI = 



-(2« + l)(ffi + M&)/2. 



(51) 
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With the help of Eq. (09 



this is readily solved to lead 
„ = [1 + §(2s + l)^ + 

Mgz)l]- X . From ijlT]). the 2k F susceptibility takes the 
form 



the simple pole expression z^f 



l + U2 S + l)(g 1 + M~g 3 )£- 



(52) 



The expected simple pole divergence at £o then occurs 
in the site CDW (A/ = +) response for the MC model 
and in the BOW(A/ = — ) response for the SSH model. 
No enhancement is found for the susceptibilities in the 
superconducting channel. 

Strictly speaking, the above adiabatic MF results hold 
for models where only momentum independent couplings 
are retained. In the case of the SSH model for spinlcss 
fermions, however, the adiabatic limit of Eqs. (|3H33[) 
does not coincide with the MF result due to the pres- 
ence of g u - In the adiabatic limit the flow equations read 



„„ 1„ 2 1,_ 

Oi9i = ~-^9\- 2 \9u + .'/;; 



7, s \2 




LO 



o 



10' 



MC 




~s=0 




f 
|g,l=0-8 


CDW — 
BOW 


1 »A= - 7 


SC 



FIG. 1: Typical variation of the susceptibilities with the scal- 
ing parameter I for the MC model for spinless fermions(s = 0, 
right) and spin-^ fermions (s = 1/2, left). The locus of the 
singularity at £ c gives the value of the gap A = Eq(£ c ). 



g u exp 



-2 / (l + ~g x {£))dt 



(53) 



The solution of these equations shows that the value of 
the adiabatic SSH gap Ao for spinless fermions is slightly 
reduced compared to the MF prediction (fSTJ)) where g u is 
absent. 



B. The molecular crystal model 

1. Spinless case 



The solution of the flow equations (|29ll30p for the MC 
model in the spinless case (s = 0) is obtained by using 
the antisymmetrized boundary conditions given in Q151 
[To]) at £ = 0. The typical flow of susceptibilities in the 
Peierls and Cooper channels at an intermediate phonon 
frequency is shown in Figure [TJ Like for the MF result 
(|52| . the singularity is found to occur solely in the site 
[M = +) CDW susceptibility at £ c . There is no no- 
ticeable enhancement of other responses including those 
of the superconducting channel. The singularity signals 
the existence of a Peierls gap A (= 2Epe~ lc ) with an 
amplitude that is reduced at non zero ujq compared to 
its adiabatic value Ao (Eq. (|5D|) ). Figure shows this 
renormalization as a function of the ratio ujq/Aq of the 
phonon frequency to the MF gap (here the molecular 
mass Mo is varied while the spring constant kq is kept 
fixed). For small loq/Aq the gap is weakly renormal- 
ized and remains close to its classical value. However, 
when the ratio ljq/Aq approaches unity, the gap under- 
goes a rapid decrease due to quantum fluctuations. This 
results from the growing of vertex corrections and inter- 
ference between Peierls and Cooper scattering channels. 



< 




FIG. 2: The site CDW gap of the MC model for spinless 
(s = 0) and spin-i (s = 1/2) fermions as a function of the 
phonon frequency and for different couplings. Both quantities 
are normalized to the MF gap. 



These fluctuations signal a change of regime (defined at 
the point of a change of curvature for the gap profile) 
that we refer to as a quantum-classical crossover for the 
gap . 

The remaining gap tail terminates with a transition 
to a A = disordered state at a threshold frequency 
slightly above Ao- The ratio ujo/Aq at which the transi- 
tion occurs is weakly dependent on the initial gi for the 
range of coupling covered by the present RG. This result 
corroborates the old two cut-off scaling arguments for 
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the disappearance of an ordered state at ljq ~ Aor^ and 
agrees with the DMRG, 9 and Monte Carle 5 results for the 
MC model. The nature of the transition to the quantum 
gapless state is also of interest. We follow the notation 
of Ref4 and define the coupling a = ^y/\gi\u}QEp. We 
see from Figure [3] that the variation of the gap A, close 
to the critical a c at which the transition occurs, follows 
closely the Baxter formula for a Kosterlitz-Thouless (KT) 
transition 3 ^ 



A 



2E P 



(54) 



where b is positive a constant. This behaviour found in 
the weak coupling range is similar to the one obtained by 
the DMRG method and perturbative expansion in strong 
coupling^ 

For phonon frequency above the threshold, the new 
state is expected to be a Luttinger liquidi&i 5 - This is 
seen at the one-loop RG level from the existence of 
a power law behaviour of the site CDW susceptibility, 
namely X^„=o(^) ~ [Eo{£)]~~ < ■ The latter takes place 
only above some characteristic i* (Fig. [4]) that depends 
on ojq and which decreases with the strength of the cou- 
pling. Non universality is also found for the Luttinger 
liquid exponent 7 for Eq(£) <C Eq(£*). Following the one- 
dimensional theor y 38 ! 39 the exponent can be written as 
7 = 2 — 2K p , where K p is the stiffness parameter for the 
density degrees of freedom that enters in the bosoniza- 
tion scheme. Within the limitation of a weak coupling 
theory, it is therefore possible to determine the depen- 
dence of K p on interaction and phonon frequency. As 
shown in Figure 03 the one-loop RG results confirms the 
non universal character of the stiffness parameter. Going 
down on the frequency scale, K p is sizeably reduced at 
the approach of the KT transition where retardation ef- 
fects have a strong influence on the properties of the Lut- 
tinger liquid parameter. We find that K p stays above the 
minimum value of \ known for isotropic spin chain in the 
gapless domain— — - following the Wigner-Jordan trans- 
formation of spins into spinless fcrmions. On the other 
hand, K p is only weakly dependent on the couplings at 
large loq, where it tends to the non adiabatic - coupling 
independent - value K p = 1 at loq — > 00. Recall that the 
initial couplings of the MC model (Eqs. (fT5|) and ifTTjl ) 
vanish in this limit, and the system is equivalent to a non 
interacting Fermi gas. 

From the variation of the critical coupling a c with the 
phonon frequency loq, one can construct the phase dia- 
gram of Fig. [6] The phase boundary between the insu- 
lating and the metallic Luttinger liquid states is found 
to follow closely a power law dependence of the form 
a c ~ u>q, with the exponent rj « 0.7. This feature 
captured by a one-loop calculation is analogous to the 
quantum-classical boundary of the phase diagram of the 
ID XY spin-Peicrls model determined by the DMRG 
method. 8 The latter model is also characterized by a zero 
temperature KT transition as we will see for the spinless 
SSH model in Sec. ITVCl 
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FIG. 3: The site CDW (left) and BOW (right) gaps as a 
function of the critical parameter of the KT transition to a 
Luttinger liquid for the CM and SSH models in the spinless 
case. The dotted line is a least squares fit to the Baxter 
formula Eq. (f54|) . 
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FIG. 4: Typical power law divergence of the site CDW (left) 
and BOW (right) susceptibilities at i > f in the gapless 
Luttinger liquid regime. 



2. Spin-^ fermions 

The results for the MC model with spin-i fermions 
(s = 1/2) ensue from the solution of Eqs. (|24H26p and the 
computation of the susceptibilities (j47| , from ()41M2[) and 
(|44M5j) . Like spinless fermions, the singularity is found in 
the M = +, site CDW susceptibility at finite l c (Fig. [J). 
The corresponding value for the gap A is reduced with 
respect to the adiabatic mean-field result Ao in (|50|) . The 
onset of quantum fluctuations due to growing interference 
between different scattering channels is again responsible 
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FIG. 5: One-loop calculation of the density stiffness parame- 
ter K p of the MC (upper panel) and SSH (lower panel) models 
in the spinless case as a function of the initial coupling \gi\, 
and for different phonon frequencies. The continuous line in 
the lower panel corresponds to the antiadiabatic one-loop re- 
sult Kp = 1 - 



for a quantum-classical crossover when ujq appoaches Ao 
where there is change of curvature in the gap profile, but 
the gap never goes to zero. It remains finite at large 
phonon frequencies and is dependent on the bare attrac- 
tive amplitude At large frequency the singularity at 
£ c occurs essentially independently for spin [gi ({u>})] and 
charge [2g 2 ({w}) — gi({u)}), <?3({w})] combinations of cou- 
plings at zero Peierls and Copper frequency. As a func- 
tion of wo, the system then undergoes a crossover from 
a renormalized classical Peierls state towards a quantum 
but still site-CDW ordered state in which both spin and 
charge degrees of freedom are gapped due to attractive 
couplings and the relevance of umklapp processes at arbi- 
trary large but finite u>q. An ordered state is well known 
to be found at large ujo in Monte carlo simulations^ This 
quantum-classical crossover marks the onset of a decou- 
pling between spin and charge degrees of freedom, a sep- 
aration found in the Luther-Emery modeL 23 i 38 ' 40 



FIG. 6: Phase diagram of MC (left) and SSH (right) models 
for spinless fermions. The full squares are the RG results and 
the continuous lines is the power law 
coupling of the KT transition on the phonon frequency. 



ujq 4 of the critical 



It is worth noting that in the purely non adiabatic 
case where luq is strictly infinite, the initial couplings 
(|13|) are independent of frequency and satisfy the con- 
ditions gi < 0, and gi — 2<?2 = \§a\, which coincide with 
those of an attractive Hubbard model. Its exact solu- 
tion is well known to give a disordered ground state. At 
the one-loop level, the RG equations (|24j|26p at zero ex- 
ternal frequencies show indeed that g\ alone is singular, 
with a gap in the spin sector only. Umklapp processes 
are irrelevant and charge degrees of freedom remain gap- 
less, consistently with the absence of long-range order at 
u>o = oo i 23 i 24 ' 25 i 41 Working at arbitrarily large but finite 
wo introduces finite retardation effect that is sufficient 
to make initial conditions deviate from those of the at- 
tractive Hubbard model. This restores the relevance of 
umklapp term in the charge sector and in turn long-range 
order £ 



C. SSH 

1. Spinless case 

We turn now to the study of the SSH model. In the 
spinless case the presence of the non local umklapp term 
g u introduces some qualitative differences with the MC 
model for which this term is absent. Thus the solution 
of (|3HI33p , and (|41I44I47|) in the spinless case shows that 
for small wn/Ao, the BOW susceptibility (/x p= o, M = —) 
is the only singular response that leads to a gap at zero 
temperature (Fig.©. 

As one moves along the uijj / Ao axis (along which the 
mass Md varies and kd is constant), one finds again from 
Fig. [51 that for lud/Aq < 0.1, the gap is weakly renor- 
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FIG. 7: Typical variation of the susceptibilities for the SSH 
model as a function of the scaling parameter I for spinless 
fermions (s = 0, right) and spin-i fermions with spins (s = 
1/2, left). 



malized compared to its adiabatic classical value com- 
puted from (|53|) . As the ratio increases further, there 
is a strong downward renormalization for the gap which 
undergoes a quantum-classical crossover. However, at 
variance with the MC model, the ratio / Aq at which 
there is a change of curvature in the variation of the gap 
shows a stronger dependence on the amplitude of the ini- 
tial coupling - parametrized by the backward scattering 
part |<7i| of gf in (fT5)) (triangles, inset of Fig. [5]). 

At higher phonon frequency we come up against a crit- 
ical value where the gap completely vanishes and the sys- 
tem enters in a metallic state at zero temperature. This 
critical ratio is also coupling dependent (squares, inset of 
Fig.[H]). In the non adiabatic limit when (wd/Ao) -1 ~~ * 0, 
the critical coupling heads on to the one-loop limiting 
value \g"i\ = 1, which can be extracted directly from 
Eq. (f33|) in this limit. A frequency dependent thresh- 
old \gf\ for the existence of an ordered state is a direct 
consequence of the relevance of the non local umklapp 
term g u in (|33p . which differs markedly from the MC 
model and simple two-cutoff scaling arguments especially 
at large phonon frequency. The phase diagram shown in 
the inset of Fig. [5] can then be obtained for the spinless 
SSH model. A critical line for (o^/Ao) -1 vs \gf\ can 
be drawn, separating the quantum disordered state from 
BOW order. The singularity of the BOW susceptibility 
in the quantum domain is shown in Fig. [7] 

The numerical solution of the flow equations is not car- 
ried out easily in the limit of very small couplings owing 
to the large number of frequencies needed to reach the de- 
sired accuracy. Our results, obtained down to \gi\ =0.1, 
tend to show, however, that the critical line extrapo- 
lates to zero at the finite value of the ratio Aq/lud « 1.1, 
which joins the value obtained from the two-cutoff scaling 
arguments. 15 Above this value, an ordered BOW state 
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FIG. 8: The BOW gap normalized to its adiabatic value as a 
function the ratio of the phonon frequency and the adiabatic 
gap for the SSH model for spinless fermions (continuous lines, 
left to right, red: \g x \ = 0.1,0.4,0.6,0.9 < blue: \g x \ = 
1.1, 1.15, 1.2 > \gi I)) and spin-i fermions (doted lines; from 
bottom to the top: \gi\ = 0.5,0.7,0.8). Inset: phase diagram 
of the spinless SSH model. Quantum BOW order/gapless 
Luttinger liquid (LL) (squares), and the quantum-classical 
crossover for the BOW order state (triangles). 



would then be found at any finite coupling. When the ra- 
tio finally crosses the quantum-classical line, the system 
enters in a Peierls BOW state similar to the one of the 
classical adiabatic limit. This quantum-classical bound- 
ary, clearly identified in Fig. [8] as a change of regime for 
the gap, is consistent with the one found by DMRG for 
the XY spin-Pcicrls chain (following the conversion of 
spins into spinless Wigner- Jordan fermions). 8 

As one moves from the quantum massive domain to- 
wards the critical line at higher frequency, the gap col- 
lapses to zero. Following the example of what has been 
done for the MC mo del, we fol low the notation of Ref. 8 
and define a c = 5'\/l.9il' J -'o-E , F, as the critical coupling 
where the gap vanishes. We thus find that close to the 
transition, A decreases to zero according to the Baxter 
expression Eq. (|5l)l for a KT transition (Fig. [3]). This 
result which carries over the whole critical line at finite 
frequency is in accord with DMRG results obtained on 
the spin-Peierls XY 8 and XXZ 10 chains. For the latter 
model a KT transition was also found by Citro et al. re- 
using the RG method in the bosonization framework. In 
the same vein, Kuboki and Fukuyama also show by per- 
turbation theory that retardation is equivalent at large 
frequency to frustration in the spin interactions^ which 
beyond some threshold is well known to promote a KT 
transition to a dimerized state.— 

As shown by the phase diagram of Fig. [5J the criti- 
cal coupling is found to follow the power law variation 
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a c ~ tOjj, with r/ « 0.7, over a sizeable range of the 
phonon frequency (deviations are found in the limit of 
small frequency). Such a power law agrees with the one 
found in DMRG for the XY spin-Peierls chain^ and is 
similar to the one obtained in the MC case for spinlcss 
fermions. 

As regards to the nature of the gapless liquid state in 
the disordered region, the situation is qualitatively simi- 
lar to the MC model. We find the presence of a non uni- 
versal power law divergence for the BOW response func- 
tion X^ p =o(2kF) oc {Eo(£)]~~t below some characteristic 
energy scale E (i*) (Fig. 01, which indicates the presence 
of a Luttinger liquid. The weak coupling determination 
of the charge stiffness parameter K p (j = 2 — 2K P ) 38 
at the one-loop level is shown in Figure [5] By com- 
parison with the MC model, K p is smaller and shows a 
stronger variation with the strength of the coupling even 
for large phonon frequencies. This is so for in the non 
adiabatic limit, the initial conditions for both <?/({w}) 
and g u ({uj}) are non zero and a massive phase remains 
possible. In this limit, the RG results join the one-loop 
relation K p = 1 — (continuous line, lower panel 

of Fig. [5]), which is known to be a good approxima- 
tion to the exact result. 37 i 43 It is worth noticing that 
although the transition remains of infinite order for any 
path that crosses the critical line from the massive sec- 
tor to the Luttinger liquid one, the characteristics of the 
latter phase, through its exponent 7 (or its stiffness co- 
efficient K p ) is strongly dependent on retardation effect. 
With the caveat of the limited accuracy of one-loop cal- 
culations for sizeable 7, the present results would indicate 
that except for the domain of large frequency, K p pene- 
trates deeply into the 'Ising sector' where K p < 1/2 at 
the approach of the critical line. 

A similar downward renormalization of K p by retar- 
dation effects has also been found by Citro et al.^ us- 
ing the self-consistent harmonic approximation and the 
RG method in the bosonization frame for the XXZ spin 
model of the spin Peierls instability. When the spins are 
converted into fermions through a Wigner- Jordan trans- 
formation, the properties of this model are encompassed 
by the flow equations (|3HI33p following a redefinition of 
the initial conditions (|15H18p . 

2. Spin-^ fermions 

The results for the SSH model with spin-i fermions 
is obtained from the solution of Eqs. (I24M26|) using the 
initial conditions (I14[) . The computation of the suscep- 
tibilities (I47|) using (|4HI42|) and (|44H45|) shows that the 
singularity and the formation of the gap remains as ex- 
pected in the BOW channel (Fig.[7|). As a result of grow- 
ing interference between the Peierls and Cooper chan- 
nels and vertex corrections in the scattering amplitudes, 
the reduction of the BOW classical gap as a function of 
the frequency, (Fig. [5]) is less pronounced than for the 
MC model. This reduction then evolves to a quantum- 



classical crossover at ujq ~ Ao . Following the example of 
the MC model, the system remains massive for both spin 
and charge and is thus BOW ordered in the quantum 
regime. The amplitude of the gap at large frequency is 
however bigger. As a matter of fact, in the non adiabatic 
case where lod is infinite, the initial couplings (fT4|) are fre- 
quency independent but at variance with the MC model, 
they satisfy the inequalities g\ < 0, g\ < g%. These are 
compatible with the Luther-Emery conditions for a mass 
in both spin and charge sectorsj 23 ' 40 ' 44 



V. CONCLUSION 

In this work we used an extension of the RG approach 
to one-dimensional fermion gas that includes the full 
influence of retardation in the interactions induced by 
phonons. Within the inherent bounds of a weak coupling 
theory, the method has been put to the test and proved 
to be rather satisfying, providing a continuous descrip- 
tion of the gap as a function of the phonon frequency 
for electron-phonon models with either spinless or spin- i 
fermions. Generally speaking, the results brought out the 
importance of the static scale Ao of the adiabatic theory 
for the occurrence of a quantum-classical crossover for 
the gap as one cranks up the phonon frequency, confirm- 
ing in passing the old arguments of the two-cutoff scaling 
approach. The RG calculations allowed us to study the 
nature of the transition to the gapless liquid phase for 
spinless fermions. For both the MC and SSH models, 
this transition was found to be of infinite order, consis- 
tently with existing numerical results. 

The RG method for the SSH model required to take 
into account the momentum dependent umklapp term. 
The latter is responsible for the continuation of the in- 
finite order critical line to arbitrary large phonon fre- 
quency where it connects to the well known results of 
frustrated spin chains. The existence in the gapless phase 
of a sharp power law behaviour of 2kp density response 
at low energy showed that this phase can be identified 
with a Luttinger liquid. The non universal variation of 
the power law exponent with the strength of interaction 
and retardation was obtained at the one-loop level. Re- 
tardation effects induce a downward renormalization of 
the Luttinger parameter K p for both models in the disor- 
dered phase. However, this renormalization is apparently 
much stronger in the SSH case where K p goes under its 
limiting i value known for the anisotropic spin chain in 
the gapless regime. 

While this work did not dwell on the combined influ- 
ence of direct and retarded interactions on Peierls-type 
instabilities, Coulomb interaction can be actually incor- 
porated without difficulty following a mere change of the 
boundary conditions for the RG flow of scattering ampli- 
tudes. A further extension of the method that includes 
both the frequency and momentum functional depen- 
dence of the scattering amplitudes would be also worth 
while. As was shown very recently by Tarn et al.^ in the 
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context of the one-dimensional MC-Hubbard model and 
by Honerkamp et al.,^ in the two-dimensional situation, 
the difficulties inherent to such an extension proved not 
insurmountable. An RG implementation of this sort for 
interacting quasi-one-dimensional electron systems would 
be quite desirable. It would yield a more complete de- 
scription of electronic phases found in correlated materi- 
als like the organic conductors and superconductors. The 
coupling of electrons to both intramolecular and inter- 
molecular (acoustic) phonon modes are in practice both 
present in these systems and their characteristic energies 
are often close to the energy scales associated to the var- 
ious types of long-range order observed^ These systems 
then fall in the intermediate phonon frequency range con- 
sidered in this work, and for which retardation effects 



may play a role in the structure of their phase diagram. 
The impact of retarded interactions on electronic states 
found in quasi-one-dimensional conductors is currently 
under investigation. 
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